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Abstract. Neutron scattering is a powerful probe of strongly correlated systems. 
It can directly detect common phenomena such as magnetic order, and can be used 
to determine the coupling between magnetic moments through measurements of the 
__i spin- wave dispersions. In the absence of magnetic order, one can detect diffuse scat- 

"Ti tering and dynamic correlations. Neutrons are also sensitive to the arrangement of 

I atoms in a solid (crystal structure) and lattice dynamics (phonons). In this chapter, 

M we provide an introduction to neutrons and neutron sources. The neutron scatter- 

ed ing cross section is described and formulas are given for nuclear diffraction, phonon 

scattering, magnetic diffraction, and magnon scattering. As an experimental exam- 
ple, we describe measurements of antiferromagnetic order, spin dynamics, and their 
evolution in the La2-i:Baa;Cu04 family of high-temperature superconductors. 
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1.1 Introduction 



^ A common symptom of correlated-electron systems is magnetism, and neu- 

^^ tron scattering is the premiere technique for measuring magnetic correlations 

T-H in solids. With a spin angular momentum of ^U, the neutron interacts directly 

^N with the magnetization density of the solid. Elastic scattering can directly 

reveal static magnetic order; for example, neutron diffraction provided the 



\l first experimental evidence for Neel antiferromagnetism [1] . Through inelastic 

^-^ scattering one can probe dynamic spin-spin correlations; in an ordered anti- 

, ferromagnet, one can measure the precession of the spins about their average 

• • orientations, which show up as dispersing spin waves. 

, 5^ Neutrons do not couple to the charge of the electrons, but instead scatter 

S^ from atomic nuclei via the strong force. Despite the name, the small size of 

$H the nucleus compared to the electronic charge cloud of the atom results in 

^ a rather weak scattering cross section. The magnetic and nuclear scattering 

cross sections are comparable, so that neutron scattering is very sensitive to 
magnetism, in a relative sense. 
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A challenge with neutron scattering is that the combination of weak scat- 
tering cross section and limited source strength means that one needs a rela- 
tively large sample size compared with many other techniques. The value of 
the information that can be obtained by neutron scattering generally makes 
worthwhile the effort to grow large samples; nevertheless, in practice it is 
useful to take advantage of complementary information obtained from tech- 
niques such as muon spin rotation spectroscopy and nuclear magnetic reso- 
nance. The latter techniques yield less information but often provide greater 
precision. There have also been continuing developments in resonant x-ray 
scattering; nevertheless, neutron scattering will remain an essential technique 
to investigate strongly correlated systems for the foreseeable future. 

As we have space only for a concise introduction to the field, we note that 
there plenty of more extended references available. A good summary of the 
theory of neutron scattering is given by Squires [2], while a more detailed 
description is provided by Lovesey [3]. We have contributed to a technique- 
oriented book [4] and to book chapters on magnetic neutron scattering [5,6], 
and new books on the subject continue to appear. 

To illustrate some of the concepts and capabilities, we will use examples in- 
volving copper-oxide compounds, especially from the family La2_a;Ba2:Cu04, 
which includes phenomena from antiferromagnetic order to high-temperature 
superconductivity. More details on neutron scattering studies of cuprates are 
given in recent reviews [7-10]. 



1.2 Basic properties of the neutron and its interaction 
with matter 

The neutron is an elementary spin-1/2 particle, which, together with its 
charged relative, the proton, is a building block of the atomic nucleus. Ac- 
cording to the "standard model" of the elementary particles, the neutron 
and proton are fermionic hadrons, or baryons, composed of one "up" and 
two "down" quarks, and two "up" and one "down" quarks, respectively. The 
basic properties of a neutron are summarized in Table 1.1. 

Table 1.1. Basic properties of a neutron. The gyromagnetic ratio, 7n, and the g- 
factor, <;„, are defined by /i„ — ^nCTn ~ —g-niiN^n, where cr„ is the neutron's angular 
momentum, S„ = (Tn/h is the neutron's spin {Sn = 1/2), /ijv = eTi/{2mpC) = 
5.05078 X 10"^^ J/T = 5.05078 x 10"^^ erg/Gs is the nuclear magneton [11, 12]. 
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Although the neutron is electrically neutral, it has a non-zero magnetic 
moment, similar in magnitude to that of a proton (/x„ — 0.684979/ip), but 
directed opposite to the angular momentum, so that the neutron's gyromag- 
netic ratio is negative. The neutron's mass, ?n„ — 1.00866 Da (atomic mass 
units) is slightly larger than that of the proton, ?ti„ = 1.00728 Da, and that of 
the hydrogen atom mn = 1.00782 Da. Therefore, outside the nucleus the free 
neutron is unstable and undergoes /?— decay into a proton, an electron, and an 
antineutrino. Although the free neutron's lifetime is only about 15 minutes, 
this is long enough for neutron-scattering experiments. For example, a neutron 
extracted through the beam-tube in a nuclear reactor has typically reached 
thermal equilibrium with the water that cools the reactor in a number of 
collisions on its way out (such neutrons usually are called thermal neutrons). 
Assuming the water has "standard" temperature of 293 K, the neutron's most 
probable velocity would be about 2200 m/s. It would spend only a fraction of 
a second while it travels along the < 100 m beam path in the spectrometer 
to be scattered by the sample and arrive in the detector. 

Neutrons used in scattering experiments are non-relativistic. Therefore, 
the neutron's energy, £■„, is related to its velocity, t)„, wave vector, fc„ = 
rrinVn/Ti, and the (de Broglie) wavelength, A„ = 27r/A;„, through 

En = Ml ^ |!M ^ ^JL^ . (1.1) 

2 2m„ 2mn\i 

Following the notation accepted in particle physics, the neutron's energy is 
measured in millielectronvolts (meV) . The neutron's wavelength and its wave 
vector are usually measured in A (1 A = 0.1 nm — 10^^ cm) and A""'^, 
respectively. Using these units, we can rewrite the Eq. (1.1) in the following, 
practical fashion: 

81 8042 
En - 5.22704 • 10"^ • vl ^ 2.07212 • kl = — - — , (1.2) 

^n 

where E^ is in meV, u„ in m/s, fc„ in A~^, and A„ in A. 

For the sake of comparison with the notations used in other techniques 
and in theoretical calculations, we list several different ways of representing 
typical neutron energies in Table 1.2. The different energy equivalents shown 
in the Table can be used interchangeably, as a matter of convenience. 



1.3 Neutron sources 

Neutrons are especially abundant in nuclei of high atomic number, where they 
can significantly exceed the number of protons. To create a neutron beam, the 
first challenge is to extract neutrons from the nuclei. The first practical source 
was the nuclear reactor, in which neutron bombardment of ^^^U nuclei induces 
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Table 1.2. Different notations used to represent the neutron's energy, e is the 
electron charge, h is the Plank's constant, c is the velocity of light, fis = e^/2meC = 
0.92740 X 10"^® J/T is the Bohr's magneton, fcs is Boltzman's constant [11]. Also 
shown are the corresponding neutron wave vector fc„ and the deBroglie wavelength 

An. 

E„ En/e E„/h E„/{hc) £;„/(2^b) E^/ke k„ \„ 

(10-19 J) (meV) (THz) (cm-^) (T) (K) (A-^) (A) 

1.60218 1000 241.799 8065.54 8637.99 11604.5 21.968 0.2860 

0.160218 100 24.1799 806.554 863.799 1160.45 6.9469 0.9044 

0.0801088 50 12.0899 403.277 431.900 580.225 4.9122 1.27909 

0.0240326 15 3.62698 120.983 129.570 174.068 2.6905 2.3353 

0.00160218 1 0.241799 8.06554 8.63799 11.6045 0.69469 9.0445 



fission, a process that releases several neutrons per incident neutron, thus al- 
lowing for a self-sustaining chain reaction. The neutrons that are released have 
a very large energy, whereas the fission cross section is enhanced by slower neu- 
trons. The slowing of neutrons can be achieved quite effectively by scattering 
from hydrogen, especially in the form of H2O, which can also act to cool the 
reactor core. In a research reactor, where one would like to extract some of the 
neutrons, the reactor moderator can be made more transparent to neutrons 
by replacing H2O with D2O (heavy water, with D representing deuterium). 
Cylindrical thimbles poking into the water moderator provide an escape path 
for neutrons, which form the beams that supply neutron spectrometers. 

Another approach is to knock the neutrons out of heavy nuclei with high- 
energy protons from an accelerator. Again, the neutrons that can escape the 
nuclei have very high energies that must be reduced by multiple scattering 
in a moderator. In contrast to a reactor, which produces neutron beams that 
are continuous in time, the proton beam provided by an accelerator can be 
pulsed, so that a spallation source typically has pulsed beams of neutrons. 
Targets can be made of a heavy metal such as tungsten, but newer sources 
with higher power tend to use liquid mercury in order to allow adequate heat 
removal. 

A list of the major operating spallation sources in the world is given in the 
upper portion of Table 1.3. Information on the available instrumentation and 
capabilities can be obtained from the listed web sites. With a pulsed neutron 
source, each burst of neutrons is produced in a narrow time window, so that 
one can distinguish between neutrons of different velocities by their travel 
time, or "time of flight" . Using a rotating shutter, one can select incident 
neutrons of a desired energy; the energy of scattered neutrons can then be 
determined by their time of arrival at a detector. 

The spallation source SINQ at the Paul Scherrer Institut provides a con- 
tinuous, rather than pulsed, beam, so its instrumentation has more in common 



1 Neutron Scattering 



o 



a 




(U 




M 




s 


rl 


a 


,0 






M 


+j 


2 


cS 







-^ 





■7^ 


h-l 









Co 





^ 




fi 




-.^ 


."^ 


AS 


(P 




> 




3 


,1^ 







> 


y 


3 


0* 















B 


"3 


1 




,ft 

ft 


ft 
■i 


a 


r^H 


fe 




& 


a 


a 


ft 


ft 


ft 


-^ 


•+^ 


+J 


-(J 


+^ 


+j 


+j 


+J 


-(J 


+J 


-a 


<< 


-B 


j: 


-a 


< 


CO 






t:3 


CO 


^ 






a 


;z) 








-5 




s" 

z 




a 


+-3 


oT 

bO 

T3 





« 
p 




CO 




< 




T3 


3 


a; 















J 


5 


rl 


> 




,i2 


XI 




r^ 


a 


cS 


<; 




a 


J 


J 




2; 


-a 

a 



S 


"ft 
ft 


bO 
H 




Q) 





T3 





Sh 






-B 


-(J 




2 


< 




CO 


^ 


S 


3 


& 


1 





J 


Cti 


^ 


Ph 



3 

o 

CO 

a 
o 

+j 
o 

o 



O 



O 
u 

3 



•S, M ^ 

n 5^ I— I Mh 

CO J £2 ^^ 



cS 



O 

1< a 



<! - 

= -- KH PL, g 

M ^ CO 



t§ 



tt; 






> 

o 









Ph Ph a a 

-^j -^j +j +j 

+j -^j +j +j 

^ ^ ^ ^ 









3 




P* 


i 


'^H 


i 


bO 

d 


% 
J 


d 


^ 


3 


^ 




^ 













;h 


M 


^ 


n3 


3 


I 


§ 


o" 


^ 


-3 



CO 
CO 

P Q 
H bS 

»■ s 

bO ,ft 






Oj 

o 



,ft 

3 
oj 

bO 

T3 



Oj 

o 



,3 
o 

cS 
0) 

0) 

3 
o 

3 

0) 



3 

a 






>1 

bO 



0) 


Z 


a 


ft 


Sh 





^ 


J 



his 


0) 


3 




i^ 


cd 


^ 


3 


J 


E 







^ 


H 


+^ 


bO 


CO 


m 



ffi z 



^ ^ 
A fe 



3 
o 



o 
,ft 



ft 







2; 











^ 




-3 


0) 




ro 


03 




oj 


oj 
ft 


-3 
bO 


3 

CO 


oj 
1-^ 




m 
3 


3" 
8 
s 




1 




■< 


CS 





3 




oj 


CO 


H 


J 




Q 







(P 
> 






3 












's 

(P 
-3 


3 
IP 

P5 


1^ 







>1 












CJ 


•••P 


oj 


< 


(P 
bo 





M 


H 




'b 


^ 


CO 


fH 


<; 


H 


rt 




-3 




< 


Oj 




CO 


H 


CO 


1 — 1 


^ 




H 


1-^ 




Z 


< 


Z 


E^ 









■< 


« 










c; 


J 


< 


< 


z 


fc 


< 





ffi 



6 Igor A. Zaliznyak and John M. Tranquada 

with reactor facilities, which arc Hstcd in the lower portion of Table 1.3. With 
a continuous source, it is common to select the desired energy of incident neu- 
trons by Bragg diffraction from a crystal (or array of crystals). In a triple- axis 
spectrometer [4], one also uses Bragg diffraction to analyze the energy and 
momentum of neutrons scattered by a sample. Again, many of the facility web 
sites provide a wealth of information on spectrometers and capabilities. 



1.4 Neutron interactions and scattering lengths 

Many of the fundamental advantages of neutron scattering techniques arise 
from the fact that the neutron's interactions with matter are usually weak and 
are extremely well understood. Hence, neutrons afford direct experimental in- 
sight into dynamical properties of the material system of interest, unperturbed 
by the probe and essentially undistorted by the details of its interaction with 
matter. These properties contrast favorably with X-ray or charged-particle 
(electron, niuon) techniques, where the probe could significantly perturb the 
system, and the interaction matrix elements between the system and the probe 
are often very complicated and profoundly impact the physics measured in the 
experiment. 

The scattering of neutrons by an atomic system is governed by two fun- 
damental interactions. The residual strong interaction (nuclear force) gives 
rise to scattering by the atomic nuclei (nuclear scattering). The electromag- 
netic interaction of the neutron's magnetic moment with the sample's internal 
magnetic fields, mainly originating from the unpaired electrons in the atomic 
shells, gives rise to magnetic scattering [2,3,13-15]. 

Magnetic interaction of a neutron with a single atom is of relativistic origin 
and is very weak, so that magnetic neutron scattering can be treated using the 
Born approximation. The interaction potential consists of the dipole-dipole 
interaction with the magnetic moment associated with the electronic spin, 
Mse = 9sSe ~ — 2se (^s ~ —2.002319 is the Lande 5-factor), 



— {fJn ■ tlse)5(r) -5 + -g 



Vse[r) = - — (Mn • Mse)<5(r) g \ g , (1.3) 



and the interaction with the electric current associated with the electron's 
orbital motion 

Vslir)^2^^B^-^'^. (1.4) 

Here hie = rxpe is the electron's orbital angular momentum, and r — r^ — rn 
its coordinate in the neutron's rest frame. 

While the neutron's interaction with the atomic nucleus is strong — the 
nuclear force is responsible for holding together protons and neutrons in the 
nucleus — it has extremely short range, < lO"'^^ cm, comparable to the size 
of the nuclei, and is much smaller than the typical neutron's wavelength. 
Hence, to describe the neutron's interaction with the system of atomic nuclei 
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in which the typical distances are about 1 A= 10^ ^'^ cm, a highly accurate 
approximation is obtained by using a delta-function for the nuclear scattering 
length operator in the coordinate representation, 

Sw(r) = &5(r„-i^^^). (1.5) 

Here Rn is the position of the nucleus and b is the nuclear scattering length, 
which is usually treated as a phenomenological parameter [16,17] that has been 
determined experimentally and tabulated [18-20]. In the Born approximation, 
the scattering length would correspond to the neutron-nucleus interaction 
described by the Fermi's pseudo-potential [21], 

VNirn,RN)^ bS{r^-RN). (1.6) 

m„ 

In general, the bound scattering length (that is, for a nucleus fixed in space) 
is a complex quantity [2, 13], b = b' — ib", defining the total scattering cross- 
section, (Ts, and the absorption cross-section far from the nuclear resonance 
capture, da, through 

a,=4^|6f aa = ^\b"\\ (1.7) 

For the majority of natural elements 6' is close in magnitude to the character- 
istic magnetic scattering length, r^ = — (g„/2)re = —5.391 fm (1 fm = 10"^"^ 
cm and r^ ~ e? /{raeC^) is the classical electron radius). 



1.5 Cross-section measured in a neutron scattering 
experiment 

In a scattering experiment, the sample is placed in the neutron beam hav- 
ing a well-defined wave vector ki and known incident flux density <l>i{ki), 
and the detector measures the partial current, 6Jf{kf), scattered into a 
small (ideally infinitesimal) volume of the phase space, cfikf = k^fdkfdf2f — 
{mnkf/h )dEfdfif, near the wave vector fe^, as indicated in Fig. 1.1. This 
measured partial current, normalized to the appropriate phase space element 
covered by the detector, yields the scattered current density. The double dif- 
ferential scattering cross-section, which is thus measured, is then defined by 
the ratio of this scattered current density to the incident neutron fiux density, 
e.g. 

d-^ajQ^E) ^ 1 5Jf{kf) 

dEdQ ^^{ki) dEdfi ' ' ' 

For each incident neutron in the plane wave state e*'"'''", the incident flux 
density is 'Pi{ki) — hki/nin- The scattered current density is determined by 
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a) 



b) 



c) 



E,k. sample 




g J. sample 




g J. sample 

i' i " 



8 = 263 
dQf 



q = kj-k^ Ej,<E^,kj, 




Fig. 1.1. Schematics of the scattering process in a neutron scattering experiment, 
(a) elastic, (b), inelastic, neutron energy loss, (c), inelastic, neutron energy gain. 



the transition rate Fi^f from the initial state \ki,S^^,rii), where the neutron 
is in the plane wave state e''"''"" with the spin S^ ^ and the scattering system is 
described by the set of variables rji, to the final state, |fe/, 5'^j, 77/). According 
to scattering theory [5,13,22], the transition rate is determined by the matrix 
elements of the transition operator (or T-matrix) T, satisfying certain operator 
equations, which depend on the scattering system's Hamiltonian, H, and its 
interaction with the neutron, V, 



r.^f 



2n 



kf,S'„j,r,f\T\KS:,,^,,n,) 



\ 2m„ 2m„ 



(1.9) 



Here E = Ef{rjf) — Ei{r]i) is the scattering system's energy gain. It is con- 
venient to introduce the scattering length operator, 6, which conveniently 
absorbs several factors, 



6(r„,S'„,ry) 



27rr 



(fe/, D„ f |r|fei, D„ 



and its Fourier transform, 6(q), 



(1.10) 



(1.11) 



Summing over all possible final scattering states, we obtain the double differ- 
ential scattering cross-section for a given initial state, \ki, S^ iiVi) i 

'^'^r&r = T: E kviH-QM)^ HEf{vf)-E,{v,)-E) , (1.12) 



■s; f.v! 



where the dependence on the spin-state of the neutron is implicit in b{—Q). 
The energy and momentum transfer to the sample are governed by the con- 
servation laws. 



Q^k,-kf, E = Ef{i^f)-E,{T^,) 



2m, 



<kf~k}). 



(1.13) 
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Finally, following Van Hove [23], one can use the integral representation of 
the delta-function expressing the energy conservation in Eq. (1.12), and the 
time-dependent scattering length operator whose evolution is governed by the 
system's Hamiltonian, 

6(g, t) = e'"*'H{q)e-'"'/^ , (1.14) 

to recast the double differential scattering cross-section in the most useful 
form of the two-time correlation function, 

"^'^ % E /^ e— *(r;,|St(-Q)6(-g, t)h,)^ . (1.15) 



dEdQ k 



S.J 



Here the sum is over all possible final spin states of the scattered neutron, S^ p 

since in the general case the scattering length operator, b{—Q,t), depends on 
the neutron spin, S'„. The sum over the final states of the sample has been 
absorbed into the expectation value of the two-time correlation function of the 
scattering length operator. The minus sign in front of Q in Eq. (1.15) follows 
from the convention adopted in the conservation laws in Eq. (1.13), where HQ 
is the momentum transfer to the sample, which is the opposite of the change 
in the neutron's momentum. The total measured scattering cross-section is 
obtained by taking the proper thermal average of Eq. (1.15) over all possible 
initial states, \r]i). 

While the scattered neutron's wave vector kf is uniquely determined by 
ki and Q, by virtue of the conservation laws (1.13), the neutron's spin state 
can be changed by transferring the angular momentum h{AS^) = ±h to 
the sample. In a polarized neutron experiment scattering between different 
neutron spin states can be measured. In such a case, the scattering length 
operator in Eq. (1.15) is a matrix with respect to different initial and final 
spin state indices; it determines the various spin-flip and non-spin-fiip cross- 
sections [3,24]. In the more common case of unpolarized neutron scattering, 
neutron spin indices should be traced out in Eq. (1.15), so that it determines 
a single unpolarized neutron scattering cross-section. 

Finally, we should mention that the double differential cross-sections in 
Eq. (1.12), (1.15) are general expressions obtained from scattering theory and 
are valid for scattering of any probe particles. The remarkable advantage of 
neutron scattering is in the fact that scattering length operators are rather 
simple, very well understood, and are directly related to the fundamental 
physical properties of the scattering sample. 



1.6 Nuclear scattering in condensed matter 

For scattering from an individual nucleus, the scattering length operator can 
be very accurately approximated by a delta-function, Eq. (1.5). For a col- 
lection of nuclei in a condensed matter system, the total scattering length 
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Fig. 1.2. Phonon-roton dispersion of the elementary excitations in the superfluid 
''He. The points show the compilation of the experimental neutron data presented in 
Ref. [25] . The solid line is the fit of the low-g part of the spectrum to the Bogolyubov 
quasiparticle dispersion. 



operator is obtained by adding scattering lengths of all nuclei, 

bN{rn) = ^ b^ S{rn - r^) , 



(1.16) 



where j indexes the nucleus at position rj with scattering length bj. For a 
system of identical nuclei, this is just a particle number density operator in 
the scattering system, times b, 



bN{rn) 



^E 

j 



H^n-rj) ^bn{rn), bN{q)^bn^ 



(1.17) 



Substituting this into Eq. (1.15) and summing out the neutron's spin states 
we obtain. 



(Pa 
dEdn 



'M 



G-"^\ih\nQn-Q{t)\'qi 



dt 
2Trh 



(1.18) 



Therefore, the nuclear cross-section measures the space-time correlation of 
the atom number density in a condensed matter system. This is exactly the 
quantity of interest in many theories of strongly-correlated quantum systems. 
One of the first successes of neutron scattering was the measurement of the 
phonon-roton dispersion of the elementary excitations in superfluid helium-4. 
Neutron data have confirmed that the shape of the dispersion is in agreement 
with that previously postulated by Landau and Feynman, as illustrated in 
Fig. 1.2. This led to the broad acceptance of the neutron scattering technique 
as a prime tool for studying quantum systems. 
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Next we consider the case in which two or more types of nuclear scatter- 
ers (with distinct scattering lengths bj and frequency of occurrence Cj) are 
present in the sample in a random fashion. For example, an element may have 
multiple isotopes, each with a distinct bj, or the nuclei have a spin, so that 
bj depends on the nuclear and neutron spin orientations, or we have at least 
two elements that are randomly distributed among equivalent positions. The 
average product of the scattering lengths for any two sites can then be written 
as 

JbJbjr) = {bril-5,,,)+¥S,,,, (1.19) 

where 

j 

^ = E^A- (1-20) 

j 

We can then distinguish between coherent scattering, 

which probes the inter-nuclear correlation, and the incoherent scattering. 



_ J^/ 

dEdfi k; 



)/'00 7, 

J2 c-^^\e-^Q-^eQ--^^*))^, (1.22) 



which probes the local autocorrelation of the nuclear position; the angle brack- 
ets denote the average over the sample state. In Eqs. (1-21) and (1.22) we have 
switched to the co-ordinate representation of nuclear density operator (1.17) 
and performed the Fourier integration. As a result, nuclear positions rj and 
r j/ (t) are quantum-mechanical operators and have to be treated appropriately 
in calculating the cross-section [2,3,14]. 



1.7 Nuclear scattering in a crystal: the Bragg peaks and 
the phonons 

In a crystal, the equilibrium positions of atomic nuclei are arranged on the 
sites of a lattice, so that the position of each individual nucleus j can be 
represented as 

r,^Rj+u,, (1.23) 

where Rj is the lattice site position, and Uj is a small displacement of the 
atomic nucleus from its equilibrium position at Rj. 

Substituting this into Eq. (1.21), one can show that the coherent nuclear 
cross-section of a monoatoniic crystal is given by 
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(1.24) 
Here ((Q • mq)^) is the time- or lattice-averaged square of the atomic displace- 
ment from equilibrium in the direction of Q, and we have taken advantage of 
the fact that the correlation function in Eq. (1.21) depends only on relative 
coordinates, which allows one summation over the N lattice sites to be com- 
pleted. The integral contains an exponentiated correlation function of atomic 
displacements. It is useful to consider the series expansion of this term in 
powers of pair displacement correlations. 

In zeroth order, the exponential factor is just 1, and one obtains the ex- 
pression for the elastic Bragg scattering in a crystal, 

i 

where we used the conventional notation for the Debye- Waller factor, with 
W = ^((Q • mq)^)- Using the lattice Fourier representation, this can be recast 
in the common form 

^ = NV* (bfc-^^ Y. ^(Q - ^) ^(^-) ' (1-26) 

where V* = {2tt)^ /Vq is the reciprocal unit cell's volume (Vb is the volume of 
the unit cell in real space), and r are the vectors of the reciprocal lattice. In 
a non-Bravais crystal, where the unit cell contains several atoms, the sum in 
Eq. (1.24) has to be split into the intra-unit cell and the inter-unit cell parts, 
leading to 

^=W*|F^(Q)p^J(Q_r)5(M, (1-27) 

where the intra-unit cell summation yields the nuclear unit cell structure 
factor, 

FNiQ) = Y. e^'^^&Me"*'^''- , (1.28) 

and fi indexes atoms in the unit cell. For some reciprocal lattice points, Fn{t) 
can be zero, which gives the Bragg peak extinction rules in a non-Bravais 
crystal. 

Expanding the exponent in Eq. (1.24) to the first order, we obtain a con- 
tribution to the cross-section that is proportional to the correlation of dis- 
placements at two different sites. In cases where static disorder is present in 
the crystal, such as dislocations or lattice strain, the time-independent cor- 
relations of displacements between different sites give rise to elastic diffuse 
scattering. 
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Calculation of the tinie-dependcnt displacements of atomic nuclei from 
their equilibrium positions in the lattice is achieved by quantizing their vibra- 
tions in terms of quantum oscillators, called phonons. A phonon is a normal 
mode of atomic vibration, a coherent wave of atomic displacements in the 
crystal. We distinguish phonons with index s. The polarization vector Eg (di- 
rection of atomic displacements) and the dependence of the energy on the 
wave vector, hujqs (dispersion), are determined by the the local inter-atomic 
potentials. The total number of such modes depends on the number of atoms 
in the unit cell of the crystal. Only three phonons, which are all acoustic, are 
present for the Bravais lattice, two transverse and one longitudinal. Taking the 
proper thermal average over the sample's equilibrium state, the contribution 
to the neutron scattering is given by 



-(S)^e-^E 



dEdn ki ^-^ 2Mujqs 

V* Y^ [5{Q q T)S{huj - hujqs){n{uj) + 1) 

T 

+ S{Q + q- T)5{huj + hiUqs)n{Lo)] , (1.29) 

where M is the mass of each nucleus. The thermal factor 

n{Lj) = {e'^^^/'^^T - l)-i (1.30) 

is the Bose distribution function describing thermal population of the oscilla- 
tor states for temperature T of the sample. The first term arises from phonon 
creation and corresponds to the neutron energy loss, while the second term 
is from an annihilation of a phonon that has been thermally excited in the 
crystal and results in the neutron energy gain. 

For a non-Bravais crystal lattice, there are also optic phonons, arising 
from the different intra-unit-cell vibrations. The total number of phonons is 
equal to Si', the number of vibrational degrees of freedom of the v atoms 
comprising the basis of the unit cell of the lattice. The contribution of each 
of these phonons to the neutron scattering cross-section is 



dEdf2 ~ ki 



6„e-^" 



X 



y* Yl ['^('3 'i T)S{huj - hujq,){n{Lo) + 1) 

T 

+ S{Q + q- T)5{hLu + hLUqs)n{iu)] , (1.31) 

where e^^ is the polarization vector for site /x in mode s. For an acoustic 
phonon in the hydrodynamic, long-wavelength (small q) and low-energy limit, 
this reduces to Eq. (1.29), where the total mass of all atoms in the unit cell, 
M — X]u -^A" should be used and one must multiply by the square of the 
elastic structure factor, \Fn{Q)\'^. 
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1.8 Magnetic scattering in a crystal. Magnetic form 
factor and spin correlations. 

The magnetic interaction of a neutron with a single atom is very weak, so 
the Born approximation provides an extremely accurate account for magnetic 
neutron scattering by the atomic electrons. In this approximation, the tran- 
sition matrix is given simply by the interaction potential, T = V , where we 
have to combine the neutron's interaction with the electron's spin and orbital 
magnetic moment, Eqs. (1.3) and (1.4). Accurate accounting for the orbital 
contribution to magnetic scattering presents, in general, a rather difficult and 
cumbersome task [3]. There are some important cases where the orbital con- 
tribution is not significant, such as transition-metal atoms in a crystal, where 
the local crystal electric field typically quenches the orbital angular momen- 
tum, or the case of s-electrons, where I = 0. Nevertheless, under some very 
general assumptions, the neutron's interaction with the electron orbital cur- 
rents can be recast in the same way as its interaction with the spin magnetic 
moment, yielding for the total magnetic scattering length. 



bAr) = -^{vUr) + Vi,{r)) 



2Trn 



Vx 



Vx'^ 



(1.32) 



where jJ^pXt) = Ms.e + l^i,e is the sum of the spin and the orbital magnetiza- 
tion associated with each electron, e. The Fourier transform of the magnetic 
scattering length (1.32), which determines the scattering cross-section, is 

^ 777 4-7r 

^™(Q) ^ -TrhTyi iP" • [Q X [g X m(Q)]]) . (1.33) 

Here m{Q) is the Fourier transform of the total magnetization density of the 
atom, 

m{Q) = ms{Q) + mL{Q) 

= / c"''^'' E (-2MB*e(5(r - r,) + /x,,^) d^r, 

•' e 

Sg is the spin operator of e— th electron, /x;.e its orbital magnetic moment 
operator. 

The cross product in Eq. (1.33) ensures the important property that only 
magnetization perpendicular to the wave vector transfer, Q, contributes to the 
magnetic neutron scattering. Adding the contributions from all atoms in the 
crystal and averaging over the neutron polarizations, we obtain the magnetic 
neutron scattering cross-section measured in an experiment with unpolarized 
neutrons {a, j3 = x,y,z), 
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(1.34) 



aJ3 



Here 

is the Fourier transformed magnetization density operator in the crystal. 
Hence, magnetic neutron scattering measures the time- and space-dependent 
correlations of the magnetization fluctuations in the sample. Introducing the 
dynamic correlation function, 

/°° rlf 

j-^-^M^M^Q{t))^-^, (1.35) 

we can rewrite Eq. (1.34) as 

where r^ = — 2/iB/i„(2?7i„/?i ) = —5.391 • 10^^'^ cm is the characteristic 
magnetic scattering length. 



1.8.1 The detailed balance constraint and the FDT 

The dynamic correlation function defined above by Eq. (1.35) obeys two im- 
portant relations that are derived in the linear response theory [2,3,15]. First, 
it is the detailed balance constraint, which relates the energy gain and the 
energy loss scattering at a temperature T, 

5"'3(Q,cj)=e^"/'=«^5^"(-Q,-w). (1.37) 

The second is the fluctuation-dissipation theorem (FDT), which relates the 
scattering intensity with the imaginary part of the dynamic magnetic suscep- 
tibility, 

Cp{Q,u) = TT (l - e-'"-/'=BT^ S"^{Q,uj) . (1.38) 

Here Xa/siQ^'^) ^-^id S"-^{Qi^) denote x^^(Q,a;) and S°'^{Q,uj) symmetrized 
with respect to {a, /3, Q} -> {/?, a, —Q}- A system with a center of inversion 
has symmetry with respect to {Q} — > {~Q}, in which case the tildes can 
be dropped for the diagonal components in {a,f3} indices. This is the case 
for which the FDT is most frequently written [5]. The FDT, Eq. (1.38), is a 
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consequence of the detailed balance condition (1.37) and the causality rela- 
tions, which require that Xq«(Q, w) is properly asymmetric. The fundamental 
laws of nature expressed in Eqs. (1.37) and (1.38) are extremely useful in 
performing and analyzing neutron scattering experiments. 

1.8.2 Elastic and inelastic scattering 

If there exists a non-zero equilibrium magnetization in the sample, (Mq) = 
{MQ(t)), where the bar over Mqlt) denotes the time-averaging, we can intro- 
duce magnetization fluctuation around this equilibrium, mQ{t) — MQ{t) — 
(Mq), and write 

S-^{Q,u;) = (M^)(MfQ)<5(M + Cci(Q,^)> (1-39) 

where the inelastic component S"Jli{Q,uj) is defined similarly to Eq. (1.35), 
but with Mq replaced by tuq. The first term here leads to elastic scattering 
which results from static magnetization in the sample, while the second term 
describes the inelastic magnetic scattering arising from its motion. Substitut- 
ing the first term into Eq. (1.36) we obtain the unpolarized magnetic elastic 
cross-section, 

where Mq is the Fourier transform of the magnetization component perpen- 
dicular to the wave vector transfer, Q. 

1.8.3 Magnetic order and magnetic Bragg peaks 

Eq. (1.40) applies equally well to all cases where static magnetism is present 
in a crystal, whether it is a long-range magnetic order leading to Bragg peaks, 
or a short-range, e. g. nano-scale magnetic correlation, resulting in an appear- 
ance of a broad magnetic diffuse scattering. In the case of a long-range order, 
the magnetization density in a crystal typically has an equilibrium static com- 
ponent, which is modulated with a wave vector Qm, 

{M{r)) =mo(r)+m(r)e*^'""'+Tn*(r)e~*^'-"', (1.41) 

where mo(r) is a real vector function that describes the ferromagnetic com- 
ponent, if present, while m{r) can be complex and describes the staggered 
magnetization. These "Bloch amplitudes" are periodic in the crystal lattice, 
and therefore can be expanded in the Fourier series, 

m(r) = -^ Vm^c*'"'', m^ = [ m(r) e"*^'^ , (1.42) 

^0 ^ Jvo 

where the integral is over the unit cell of the nuclear (paramagnetic) crystal 
lattice. 
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Substituting (1.41) and (1.42) into Eq. (1.40), we obtain the following 
expression for magnetic Bragg scattering associated with the long-range mag- 
netic order at a wave vector Qm, 



"^'""^'^ ^ NrlV* 



dEdQ 



mt 



2Mi 




[5{Q - Q 



6{huj). (1.43) 



Here the summation is over the paramagnetic crystal lattice. This is the "large 
Brillouin zone" description, which is the most general one, in that it does not 
rely on the existence of a commensurate magnetic superlattice with a unit 
cell containing some integer number of nuclear lattice unit cells, and applies 
to incommensurate, as well as commensurate magnetic structures. Such a 
description is most convenient for stripe phases in the cuprates, which are 
often incommensurate. 

The intensities of magnetic satellites, |97i:^p, are given by the Fourier am- 
plitudes of the magnetization, (1.42), which are obtained by performing the 
Fourier integrals over the unit cell of the paramagnetic lattice. In the case 
where the unit cell magnetization could be approximated by a number of 
point-like magnetic dipoles fii, located at positions r,y, these amplitudes be- 
come the conventional unit cell magnetic structure factors, 



in{r) 



i.i' 



H„S{r - Rj - r^) , 



^^^i'^'^ 



(1.44) 



In discussing magnetic scattering we assume a rigid lattice, neglecting atomic 
displacements due to disorder and vibrations discussed above. The leading 
correction to this description is obtained by multiplying expressions for mag- 
netic cross-section with the Debye- Waller factor, e" 



-2W 



1.8.4 Magnetic form factor and spin correlations 

In many important cases the magnetization density in the crystal is carried 
by electrons localized on atomic-like orbitals, which are specified by the local 
atomic variables, such as spin and orbital quantum numbers. In such cases, 
the matrix element of the atomic magnetization in the magnetic neutron scat- 
tering cross-section can be factorized into the product of the reduced matrix 
element (form factor), which does not depend on the direction of the atom's 
angular momentum quantum numbers, and the Wigner 3j-symbol, which en- 
tirely accounts for such dependence. Hence, the cross-section can be expressed 
in terms of a product of the Q— dependent form factor, which accounts for the 
shape of the magnetization cloud associated with the atomic spin and orbital 
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variables, and a dynamical correlation function between these local angular 
momentum variables at different lattice sites. 

For magnetic ions obeying Hund's rule, neutron scattering usually probes 
states belonging to the same multiplets of the angular momentum, AL = 0, 
AS = for the Russel-Saunders atoms with weak spin-orbit and strong crystal 
field, or AJ = for the case of strong spin-orbit coupling, where the total 
angular momentum J ~ L + S is a good quantum number, such as in rare 
earths. Hence, we can write for the Fourier transform of atomic magnetization, 

{Vf\M{Q)\7j,) = -2fiBFs{Q){vf\S\m) - ^iBFL{Q){vf\L\v^) , (1.45) 

where the spin and the orbital magnetic form factors are, 

(T]'f,L, S\ V^ c-"^"-' (se • S) Iry', L, S) 
FsiQ) = ^^ '^^ ^(5+1) ' (l-^^) 

Mf, L, S\ V^ e-"^"-" (fie I ■ L) M,L, S) 

Fl{Q) - ^'^ '^^ ,,,^' '^'^ ' , (1.47) 

fiBL(L + 1) 

where we made explicit that initial and final states of the sample belong to 
the same L and S multiplet. Similar relations hold for the J multiplet in the 
strong spin-orbit coupling limit. 

Typically it is possible to define an effective spin operator, 

{7jf\M{Q)\rj,) = -g^^BF{Q){vf\S\v^) , (1.48) 



Wf.,L,S\iM{Q)-S\\vi,L,S) 

F{Q) = — ^ ~,~ , = -FsiQ) + ^^Fl(Q) . 

g^iBS{S+l) g 9 

(1.49) 

where g and gs are the effective 5— factors, (77/li + 2S'|?7i) — g{rif\S\rii), 

(77^|2S'|?7i) = gs{Tif\S\rii). These expressions are exact in the cases of a J- 

multiplet, where S ~ J, gs = '2{g — 1) and (/l = 2 — g, or a pure spin 

multiplet, where g ~ gs and the orbital contribution is absent. They give the 

leading-order approximation in other cases. If the orbital moment is nearly 

quenched, as it is for magnetic d-elements in strong crystal field, then 5' « S, 

gs ~ 2, and the orbital contribution to F{Q), is small. Assuming this to be 

the case, we shall omit tildes and use S for the effective spin. 

Using the factorization of atomic magnetization provided by Eqs. (1-48) 

and (1.49), the magnetic neutron scattering cross-section (1.34) can be recast 

as 



(1.50) 



d'^arn 


h-2X-f, QcQA^^ F*{Q) F,,{Q) 

/•OO I1 

/ e-'"*e-*Q-(«^-«.')(5"5f,(i))^ ' 
J-00 27r?i 


dEdQ 
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where we allow for the possibility that the g— factor is anisotropic, and that 
both gaj and Fj{Q) could be different for different sites j,j' of the lattice. 
Equation (1.50) relates the magnetic cross-section to the dynamic spin struc- 
ture factor, which is the Fourier transform of the time-dependent two-point 
correlation function of the atomic spin variables on the sites of the lattice, 

S'^^{Q, Lj) is a quantity which is calculated in theoretical models based on the 
local spin Hamiltonians. It also obeys a number of important relations, known 
as sum rules, which are extremely useful in analyzing neutron scattering data. 
The zero moment sum rule is obtained by integrating Eq. (1.51) in Q and cj, 
providing the direct connection of the integral neutron intensity with the spin 
value S in the lattice spin Hamiltonian, 



E 



/oo 
S°'"{Q,oj)d^qd{huj)^S{S + l). (1.52) 

-oo 



The first moment sum rule relates X^q, /_ TiujS°'°'{Q,uj)(Pqd{huj), which is 
the integral oscillator strength of the fluctuation spectrum, with the bond 
energies in the spin Hamiltonian, and so on. 

1.8.5 Spin waves 

Representing the neutron scattering cross-section via two-point dynamical 
spin correlation function, as in Eq. (1.50), is possible in a large number of im- 
portant magnetic systems, such as cuprates and other 3d magnetic insulators. 
Such a representation is extremely useful, as it allows one to connect the mea- 
sured magnetic neutron intensity with the theoretically predicted properties 
of model spin Hamiltonians, such as the Heisenberg spin Hamiltonian, 

J J' Q 

Here Jjj' = Ji^jj') is the exchange coupling between sites j and /, and Jq 
and Sq arc the lattice Fourier transforms, 

■h = T. Jn'^'"'""' ^ Sq = J2 S.o-'"-' . (1.54) 

r ■ f i 

3J -I 

In many systems with magnetic order, the average value of spin at each 
lattice site in the ground state (GS) is "frozen" at nearly the full saturation 
value, {S^) « S. In particular, this is a very good approximation for the semi- 
classical spins, S* ^ 1, in more than one dimension (ID). For quantum spins, 
S = 1/2, and/or in the low-dimensional, or frustrated systems, the order may 
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be weak, or absent, and such a picture is inadequate. Nevertheless, in a large 
number of systems magnetic order in the ground state is well developed, and 
the semiclassical spin-wave picture applies. 

Spin excitations in a magnetic system with a well-ordered ground state, 
such as a ferromagnet, where all spins are parallel, or a semi-classical anti- 
ferromagnet, where there are two antiparallel sublattices, can be visualized 
as small oscillations of classical spin vectors around their equilibrium posi- 
tions in the GS spin structure. Their wave-like spatial composition results 
from the translational symmetry of the system. Frequencies of such spin- wave 
oscillations can be calculated from the spin Hamiltonian, such as (1.53), to- 
tally within the classical mechanics, simply by writing the torque equations 
of motion for the classical spin angular momenta. For example, in the case 
of the Hcisenberg Hamiltonian (1.53) for a magnetically ordered system char- 
acterized by the ordering wave vector Qo (this includes ferromagnetism with 
Qo = 0, as well as antiferromagnetism and helimagnetism) , one obtains the 
spin- wave dispersion [26,27], 



?.c.,^2^./(J,-jQj(' -^'^+Q"y^-Q" -jQj. (1.55) 



This can be recast as Wq — ^JlJquJq^, where Tiloq = 2S{Jq — Jqo)- 

Spin waves are the normal modes of the linearized equations of motion. 
They involve small spin deviations that are perpendicular to the equilibrium 
spin direction. Hence, spin waves are transversely polarized, with two mutually 
orthogonal linear polarizations of spin oscillations possible. For a spin system 
on a Bravais lattice there are two spin-wave modes. 

In a quantum-mechanical treatment of spins, the spin-wave calculation 
proceeds via an approximate mapping of spin operators to Bose creation- 
annihilation operators, i.e. to local oscillator modes. Hence, the so obtained 
spin-wave theory (SWT) describes spin excitations as coherent waves of small 
oscillations around the local equilibrium positions, in many ways similar to 
phonons. The resulting expression for the spin-wave contribution to the neu- 
tron magnetic scattering cross-section in a sample with a spiral spin structure 
with the propagation vector Qg is 
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(1.56) 



where z is the direction normal to the plane of the spiral, and we have re- 
stricted consideration to the case of a Bravais lattice and retained only the 
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contribution corresponding to creation of a single spin wave. The contribu- 
tion arising from the absorption of a spin wave is written similarly to that 
of a phonon in Eq. (1.29). For a ferromagnet, the single spin-wave magnetic 
cross-section simplifies to, [2,3], 

dEdQ fc/" l2^^^ 

V*^\5{Q-q- T)S{huj - hujg){n{uj) + 1) + 

T 

S{Q + q-T)S{huj + hujq)n{uj) , (1.57) 

where Qu is the wave vector component along the ferromagnetic ordered mo- 
ment and we have retained the contributions from both the creation and the 
absorption of a spin wave. 

1.8.6 Anisotropic magnetic form factor and covalency 

It is clear from Eq. (1.50) that even the exact knowledge of the dynamical 
spin structure factor (available from theory in some special cases, such as in 
one dimension) is insufficient to reproduce the measured magnetic scattering 
cross-section. One also has to know the magnetic form factor, which needs 
to be obtained from an ab initio calculation of the electronic density in the 
crystal. 

In the most common case of a Hund's ion with 25" unpaired electrons 
forming spin (25'-|-l)-multiplet, the spin magnetic form factor (1.46) becomes 

2S 2S 

FsiQ) = ^Y.fc-^'^-'^\Mr)\'d'r=^J2^sAQ), (1-58) 

e— 1 e— 1 

where the sum is only over the unpaired electrons. The single-electron density, 
|'!/'e(''')p7 is determined from the many-electron atomic wave function through 
IV'e(''')P = {il'y L, S\S{r — Te)|?7', L, S). The magnetic form factor for an atom 
is therefore simply an average of those for each of the unpaired electrons. 
Similarly, the orbital form factor is the Fourier-transformed average density 
of the uncompensated orbital currents in the atom. 

If the average Hartree-Fock potential acting on an unpaired electron e 
in the atom is spherically symmetric, then the effective one-electron wave 
functions in (1.58) are the eigenfunctions of angular momentum and are 
tagged by the n,l,m = P quantum numbers, V'e('") = V'n,/,m('')- The angular 
and the radial dependencies of the electronic density factorize, 4'n,i.m{f') = 
Rn.i{'r)Yf^{6^(t)), where Y[^{6,(f)) is the spherical function giving the depen- 
dence on the polar angles 9, 4>. This so-called central field approximation is 
good when the contribution to the potential from electrons in the incomplete 
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shell is small. However, it also becomes exact for an almost-filled shell, with 
only a single electron, or a single hole, as in the case of Cu^+, or for a nearly 
half-filled shell, because the average potential of the closed, or half-filled shell, 
is spherically symmetric. 

In the general case, a single-electron wave function can always be expanded 
in a series in spherical harmonics. In each term of such an expansion, the ra- 
dial and the angular parts are again factorized, and the magnetic form factor 
is a sum of Fourier-transformed terms with different I and m. The same kind 
of an expansion is encountered in calculating the orbital contribution to the 
magnetic form factor. This is known as a multipole expansion [3]. The calcu- 
lations are ion-specific and extremely cumbersome. The general expressions 
can be obtained only for the leading, isotropic contributions, in the limit of 
small wave vector transfer, known as the dipole approximation, 

FsiQ) - UoiQ)) , Fl{Q) = \ (Oo(Q)) + {n{Q))) , (1-59) 

where jo{Q) and J2{Q) are the l,m dependent radial integrals quantifying 
the radial wave function [5,18]. The radial integrals for most known mag- 
netic atoms and ions have been calculated numerically from the appropriate 
Hartree-Fock or Fock-Dirac wave functions and are tabulated in Ref [18]. The 
fuU F{Q) is given by Eq. (1.49). 

Although the dipole approximation (1.59) is the one most commonly used, 
it is extremely crude. In particular, it does not account for the anisotropy of 
the magnetic form factors, which can be very important for ions with only one 
or two unpaired electrons. The anisotropic magnetic form factor of a single 
5d hole in a t2g orbital of the magnetic Ir^^ ion in the cubic K2lrCl6 was 
studied in Ref. [28] . The authors found that the anisotropy of the magnetic 
form factor is very large, with an additional enhancement coming from the 
hybridization of the Ir 5(i-orbital with the CI p-orbitals. 

The anisotropy of the magnetic form factor is also very pronounced in 
La2Cu04, YBa2Cu30e+y, and related cuprate materials, including the high- 
Tc superconductors, where in the ionic picture a single unpaired magnetic elec- 
tron occupies a idx^-y^ orbital. In Ref. [29] the authors found that properly 
accounting for the anisotropy of the Cu^^ magnetic form factor is essential for 
understanding the magnetic Bragg intensities measured in YBa2Cu306+y at 
large wave vectors, and can also explain the peculiar Q-dependence of the in- 
elastic magnetic cross-section in this material. Accounting for the anisotropic 
Cu^+ form factor was also very important in analyzing neutron scattering by 
high-energy spin waves in La2Cu04 [30,31], and the chain cuprates SrCu02 
and Sr2Cu03 [32,33]. The magnetic excitations in these cuprate materials 
extend to several hundreds of meV. Consequently, the measurements require 
very large wave vector transfers, for which the anisotropy of the magnetic 
form factor is very pronounced. 

The ionic magnetic form factors for 3(i orbitals can be explicitly computed 
by Fourier transforming the corresponding spherical harmonics. In particular. 
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for the (1^2 _y2 orbital relevant for Cu^+ one obtains [5], 

F{Q) = (MQ)) - y(.72(g)> (1 - cos2 Oq) 

+ ^ (MQ)) (l - 10 cos2 9q + ^ cos* 9q\ 

+ y 04(Q)) sin* Oq eos (#q) , (1.60) 

where 6q, (J)q are the polar angles of the wave vector Q in the local coordinate 
system used to specify the proper orbital wave functions in the crystal field. 
Although using the anisotropic ionic magnetic form factor of Cu^+ is much 
better than using a spherical form factor of the dipolc approximation, it is still 
not sufficient for cuprates, as it neglects the effects of covalency (i.e., charge 
transfer to the neighboring oxygen) that are expected to be very significant in 
these materials. In Ref. [33] it was discovered that covalent bonding results in 
a marked modification of the magnetic form factor in the quasi-lD antiferro- 
magnet Sr2Cu03. The local structure of the planar Cu-0 square plaquettes 
in this material is essentially identical to that in La2Cu04. Making use of 
a precise theoretical result for the excitation spectrum available in ID, the 
authors demonstrated that a good fit to the data requires a form factor that 
takes account of hybridization between the half-filled Cu 'idj.2_y2 orbital and 
the ligand O 2pn orbitals, as given by a density functional calculation. The 
hybridization causes the spin density to be extended in real space, resulting 
in a more rapid fall off in reciprocal space compared to a simple Cu^+ form 
factor, as illustrated in Fig. 1.3. Smaller values of magnetic form factor at rel- 
atively large wave vectors, where the measurement is performed, lead to the 
suppression of magnetic intensity, which could be as large as a factor of two or 
more [33] . Finally, we note that a study of covalent NMR shifts by Walstedt 
and Cheong [34] found that barely 2/3 of the spin density in La2Cu04 resides 
on the copper sites, in excellent agreement with the Sr2Cu03 neutron data of 
Walters et al. [33]. 



1.9 Application to cuprate superconductors 

The discovery of high-temperature superconductivity in La2-i;Ba2;Cu04 
(LBCO) came as a considerable surprise [35], as ceramic oxides were gen- 
erally considered to be poor conductors. The structure of LBCO and related 
cuprates involves Cu02 layers, with the Cu atoms forming a square lattice 
with bridging O atoms, as shown in Fig. 1.4(a). Anderson [36] predicted that 
the parent compound, La2Cu04, should have strong antiferromagnetic (AF) 
superexchange interactions between nearest-neighbor Cu atoms. The occur- 
rence of antiferromagnetic order was demonstrated by Vaknin et al. [37] us- 
ing neutron diffraction on a powder sample of La2Cu04. As illustrated in 
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Fig. 1.3. Wave vector dependence of the ionic magnetic form factor of Cu^^ given by 
Eq. (1.60) (solid line) and the covalent magnetic form factors for Sr2Cu03 (dashed 
line) and La2Cu04 (dash-dotted line) obtained from the ab initio density functional 
calculations [33]. Panels (a) - (c) show the dependence along three principal direc- 
tions. 



Fig. 1.4(b), the antiferromagnetic Neel order doubles the size of the unit cell 
in real space, which results in magnetic superlattice peaks, as shown in (c). 
Thus, the antiferromagnetic order can be detected through the appearance 
of superlattice peaks. The challenge in this case is that one must distinguish 
from structural superlattice peaks due to staggered rotations of CuOg octa- 
hedra [38]. Fortunately, the AF and structural peaks appear at inequivalent 
positions. 

The ordered pattern of the octahedral tilts is associated with an or- 
thorhombic distortion of the crystal structure that makes the diagonal di- 
rections of a Cu-0 plaquette inequivalent [38], as indicated in Fig. 1.5. By 
analyzing the Q dependence of the AF Bragg peak intensities, it was possible 
to determine that the magnetic moments on Cu atoms lie within the Cu02 
planes, pointing along the orthorhombic b axis [37]. Furthermore, it was pos- 
sible to show that the relative arrangement in neighboring planes is as shown 
in Fig. 1.5. With the magnetic structure determined, one can evaluate the 
magnitude of the magnetic moments by normalizing the AF peak intensities 
to the nuclear intensities and correcting for the magnetic form factor. Early 
studies yielded a small ordered moment whose magnitude was correlated with 
the magnetic ordering, or Neel, temperature, Tn [39] . Neutron scattering stud- 
ies on carefully prepared single-crystal samples eventually demonstrated the 
impact of interstitial oxygen, within the La2 02 layers [40]. Removing the ex- 
cess oxygen by annealing, one can achieve Tn — 325 K [41] and a magnetic 
moment of 0.60 ± 0.05 /ie [39]. 

If one assumes a g factor of roughly 2, then the ordered moment yields 
an average ordered spin {S) ~ 0.3, compared to the expected S — ^ per 
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Fig. 1.4. (a) Structure of a Cu02 plane, with Cu atoms indicated by filled circles 
and O atoms by open diamonds, (b) Schematic of antifcrromagnctic order, with al- 
ternating up (filled circles) and down (open circles) spins. The solid line indicates the 
chemical unit cell, while the dashed line indicates the doubled area of the antiferro- 
magnetic unit cell, (c) Reciprocal space showing fundamental Bragg peak positions 
(filled circles) and antifcrromagnctic superlattice peak (open circle) at (|, |). 



Cu atom. The reduction results from the strongly anisotropic structure and 
the low value of the spin. For a two-dimensional magnetic system described 
by a Heisenberg spin Hamiltonian, long-range order is destroyed at any fi- 
nite temperature by thermal excitation of spin fluctuations. For La2Cu04, 
weak (nearly-frustrated) couplings between the planes enable the ordering at 
finite temperature [43]. Nevertheless, the spin correlations have a strongly 
two-dimensional (2D) character, as demonstrated by neutron scattering stud- 
ies [44] . The small magnitude of the spin, combined with the enhanced zero- 
point spin fluctuations in 2D, puts the system close to a quantum critical 
point [45]. Although the large fluctuations cause problems for perturbation 
theory, spin-wave theory nonetheless yields a result, (S) = 0.3, that is very 
close to the value obtained from experiment [46] . 

The exchange couplings between the spins can be determined by analyzing 
the dispersion of the spin excitations, which can be obtained by inelastic 
scattering measurements on a single-crystal sample. Early studies of La2Cu04 
with triple-axis spectrometers demonstrated that the superexchange energy J 
coupling nearest-neighbor spins is greater than 100 meV, and that effects such 
as exchange anisotropy and interlayer coupling are very small [43]. Time-of- 
flight techniques were required to measure the highest-energy spin waves [47] , 
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Fig. 1.5. Structure of La2Cu04, with arrows indicating the arrangement of the 
magnetic moments in the antiferromagnetic state. Figure reprinted with permission 
from Lee et al. [42]. Copyright (1999) by the American Physical Society. 



and these have been refined over time [30,31]. The most recent results, from 
Headings et al. [31], are shown in Fig. 1.6. The line through the data points 
corresponds to a fit with linear spin-wave theory, which works surprisingly 
well in light of the the large zero-point fluctuations. One impact of the latter 
is the renormalization factor, with a fitted value of Zd = 0.4 ± 0.04, that 
is required to fit the measured intensity. This is somewhat smaller than the 
value Zd w 0.6 that is predicted from quantum corrections to linear spin 
waves [48]. It should be noted, though, that this analysis did not take account 
of the hybridization effects on the magnetic form factor, discussed in Sec. 1.8.4, 
which would account for some of the apparent renormalization. An additional 
effect is the damping and broadening of the energy-dependence of the spin- 
wave line shape at the zone boundary position Q = (5,0). This appears to be 
the result of interactions with a multi-magnon high-energy continuum [31]. 

The fitted dispersion corresponds to J = 143 ± 2 meV, but also requires 
longer-range exchange couplings — second and third neighbor couplings J' and 
J", which are relatively weak, and a significant 4-spin cyclic exchange term 
Jc ~ 0.4 J. The overall bandwidth of the magnetic spectrum is ~ 2 J. A recent 
analysis of the couplings, including Jc, in terms of a single-band Hubbard 
model has been given by Dalla Piazza et al. [49]. 

To achieve superconductivity, one must dope charge carriers into the Cu02 
planes. Substituting Ba^+ or Sr^+ for La^+ introduces holes. A small density 
of holes, p « 2%, is enough to kill the long-range AF order, which is followed 
by a regime of spin-glass order [7]. Doping beyond p ~ 0.055 yields supercon- 
ductivity. The maximum superconducting transition temperature T^ occurs 
for p '^ 0.16, with Tc heading towards zero for p > 0.25. Inelastic neutron 
scattering studies have been performed on single crystal samples across this 
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Fig. 1.6. Spin wave (a) dispersion and (c) intensity measured in antiferromag- 
netic La2Cu04 at T = 10 K. Lines through data correspond to fits with spin-wave 
theory; the fit to the intensity includes a renormalization factor Zd = 0.4 ± 0.04. 
Figure reprinted with permission from Headings ei al. [31]. Copyright (2010) by the 
American Physical Society. 



entire doping range [7,9]. A couple of the key features are: 1) the bandwidth of 
strong spin-fluctuation scattering decreases linearly with doping, being quanti- 
tatively similar to the pseudogap energy extracted from various electron spec- 
troscopies [9,50], and 2) the wave vector characterizing the low-energy spin 
excitations splits about the AF wave vector, becoming incommensurate [7,51]. 
Insight into the cause of the magnetic incommensurability was provided by 
neutron diffraction measurements on a closely related material, 
Lai.48Ndo.4Sro.i2Cu04 [52]. The impact of the Nd substitution is to mod- 
ify the tilt pattern of the CuOe octahedra such that the in-plane Cu-0 bond 
directions become inequivalent [38]. New superlattice peaks were observed in 
this low-temperature phase, with in-plane wave vectors Q = (| ± e, i) and 
(|, I ± e) corresponding to spin order and (±2e,0) and (0,±2e) associated 
with modulations of atomic positions due to charge order, with e « 0.12. Such 
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Fig. 1.7. Cartoons of equivalent domains of (a) vertical and (b) horizontal bond- 
centered stripe order within a CUO2 plane (only Cu sites shown). Note that the 
magnetic period is twice that of the charge period. The charge density along a stripe 
is one hole for every two sites in length. The registry of the stripes with respect to the 
lattice (for example, site-centered vs. bond-centered) has not yet been determined 
experimentally. 



results have been confirmed in the system La2-£!;Baa;Cu04 [53,54]. Analysis 
of the superlattice peaks indicates that they are evidence for spin and charge 
stripe order [55,56], as illustrated in Fig. 1.7. Because of the crystal symmetry, 
the orientation of the stripes rotates 90° from one layer to the next. 

The occurrence of maximum stripe order corresponds to a strong suppres- 
sion of the bulk Tc at p « | [54,57], suggesting that stripe order competes 
with superconductivity; however, recent studies have demonstrated that 2D 
superconductivity can coexist with stripe order [58]. It now appears that su- 
perconducting order can intertwine with stripe order [59]. Thus, understand- 
ing stripe correlations may provide valuable insights into the nature of the 
superconducting mechanism of cuprates. 

Neutron scattering on a time-of-fiight instrument has been used to char- 
acterize the spin excitation spectrum in La2_2;BaxCu04 with x = 1/8 [60]. 
The effective dispersion and the Q-integrated spectral weight are shown in 
Fig. 1.8. Above 50 meV, the excitations disperse upwards like antiferromag- 
netic spin waves with an energy gap; the solid line through the points in each 
panel corresponds to a two-leg spin ladder model with J — 100 meV. Below 
50 meV, the excitations disperse downwards toward the positions of the in- 
commensurate magnetic superlattice peaks. When the sample is warmed to 
a state with no static stripe order, the spectrum maintains its essential fea- 
tures [53,61]. It appears that stripes, whether static or dynamic, provide a 
way for the superexchange mechanism to survive when the antiferromagnetic 
layers are doped with holes. 

The relevance of charge-stripe order is less clear in cuprates families such 
as YBa2Cu3064.y and Bi2Sr2CaCu208+5; nevertheless, the dispersion of the 
magnetic excitations in these compounds (measured by neutron scattering) 
has been shown to be quite similar to that of LBCO [8,9]. The main difference 
is that the low-energy excitations tend to be gapped in the superconducting 
state, with a pile up of weight ("resonance" peak) appearing above the gap 
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Fig. 1.8. (a) Q-integrated spectral weight and (b) effective magnetic dispersion in 
the stripe-ordered phase of La2-a;Ba2:Cu04 with x = 1/8, from [60]. The solid lines 
through the data points are described in the text. In (a), the peak at ~ 40 meV is 
now know to be due to a phonon mode. In (b), the effective dispersion is plotted for 
q along a line through the incommensurate magnetic superlattice peaks. 



for T < Tc- The commonality of the dispersions over a broad energy range 
suggests that the charge and spin correlations in superconducting and striped 
cuprates are similar. 
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